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Abstract 

A static extended trailing edge attached to a NACA0012 airfoil section is studied for achieving lift enhancement at a 
small drag penalty. It is indicated that the thin extended trailing edge can enhance the lift while the zero-lift drag is 
not significantly increased. Experiments and calculations are conducted to compare the aerodynamic characteristics 
of the extended trailing edge with those of Gurney flap and conventional flap. The extended trailing edge, as a 
simple mechanical device added on a wing without altering the basic configuration, has a good potential to improve 
the cruise flight efficiency. 
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Introduction 

It has long been recognized that flaps can significantly alter wing aerodynamics for high lift generation. 
Conventional flaps have larger drag; therefore they are mainly deployed for takeoff and landing and are not suitable 
in cruise flight. Small trailing edge devices like Gurney flaps and divergent trailing edges for lift enhancement have 
attracted considerable attention in aeronautical community. The Gurney flap was introduced by Liebeck [1] for 
aeronautical applications, and considerable measurements and calculations have been performed to determine the 
aerodynamic characteristics of wings with Gurney flaps at low speeds [2-9]. Measurements were also made at high 
speeds [10]. Gurney- flap-type tabs have been used along with vortex generators to enhance lift and suppress flow 
separation for a single-slotted-flap airfoil at large flap deflection [11]. Microflaps that can change the deflection 
angle were used for flutter suppression [12]. To mimic the behavior of bird feathers during landing, a self-adjusting 
movable flap close to the trailing edge was proposed by Bechert et al. [13], which was able to maintain higher lift 
when stall occurs. A review on trailing edge devices and other control technologies is given by Stanewsky [14]. 

Measurements of avian wing geometry for extracting typical avian airfoil sections indicate that the merganser and 
owl wings have a very thin trailing edge that is a single layer of feathers extended from a ‘normal’ airfoil section 
[15]. Naturally, a question is whether and how this unique trailing edge affects the global aerodynamic 
characteristics of the wing. This inspires the concept of a static extended trailing edge (SETE) as illustrated in 
Figure 1(a). The main airfoil element remains unchanged, but it is extended at the trailing edge by attaching a thin 
splitting plate of suitable length and rigidity. As shown in Fig. 1(b), the basic geometrical parameters are the main 
airfoil chord (c), length of SETE (/) and deflection angle ( S ). Depending on specific applications, the extension 
could be an aluminum plate, polymer membrane, composite sheet or smart material plate. In general, a thin 
extended trailing edge is flexible such that it can be passively changed through flow-structure interaction or actively 
controlled by embedded actuators. In the current operating conditions, our analysis shows that when extended 
trailing edges are sufficiently rigid, their shapes are not sufficiently affected by flow and no vibration is induced. 
Thus, the cases studied here are static or quasi-steady. The effects of SETE on the wing aerodynamics are mainly 
due to modifications of the airfoil camber and of the flow structure at the trailing edge. The lift enhancement by 
SETE is expectedly due to the camber effect. It is speculated that the drag penalty of SETE is small since it is 
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usually embedded in the wake of the main airfoil. The drag penalty should be examined carefully in comparison 
with Gurney flaps and conventional flaps. The objective of this work is to study SETE as an alternative flow control 
technique for lift enhancement at a small drag penalty in cruise flight. 

We first describe the geometrical features of SETE and relevant similarity parameters. Experimental results are 
presented including the lift and drag coefficients of the baseline NACA0012 airfoil model and model with SETE and 
Gurney flaps. The lift enhancement is achieved by SETE, while the drag polar and lift-to-drag ratio remains largely 
unchanged. The zero-lift drag and Oswald efficiency are not significantly affected by SETE. The benefit analysis 
for SETE in cruise flight is given in comparison with Gurney flap. Based on thin airfoil theory and CFD 
simulations, the aerodynamic characteristics of SETE are further compared with those of the corresponding Gurney 
flap and conventional flap. 



(a) (b) 

Figure 1. (a) Illustration of an extended trailing edge, (b) Geometrical parameters and coordinate system. 

The team members of experimentalists and computationalists are fully engaged throughout the various phases of the 
research. The open and continuous communications have facilitated the planning, the design, and the execution of 
experimental and computational tasks. In fact, they often lead to valuable insights that are not readily apparent to a 
uni-approach team. We believe a full integration of experimental and numerical approaches can be very effective 
and cost-efficient for aerodynamics studies. 

Wind Tunnel and Model 

Tests were conducted in the Advanced Design Wind Tunnel (ADWT) at Western Michigan University (WMU). 
The ADWT is a low speed, closed circuit single return, atmospheric tunnel with a test section of 32 inches high, 45 
inches wide, and 8 feet long. Full test section length Plexiglas windows on sides and ceiling allow optical access. 
The 161 feet long circuit consists of two legs of 65.6 feet long and the other two legs of 14.9 feet long. Low 
turbulence and good flow direction control are obtained by a combination of features, including four screens and 
honeycomb in the stilling camber, a 10:1 contraction ratio between the stilling chamber and test section, empirically 
tailored test section walls, tangential flow breathing at the test section exit to the diffuser, and a 5.4° diffuser angle. 

Under the present test conditions at Re c = 4.74 xlO 5 , natural boundary layer transition was observed at about 30%c 
on the NACA0012 airfoil model at zero AoA based on global luminescent oil film skin friction measurements [16]. 
The main drive power source is a 125 hp DC motor with a solid-state speed controller. The 16-blade fan is 6 feet in 
diameter. The flow velocity range is approximately 20 to 240 ft/sec (6 to 73 m/sec) where the total temperature is in 
the range of 60 to 120° F. This corresponds to the dynamic pressures up to 65 lb/ft 2 and Reynolds numbers up to 1.3 
x 10 6 per foot. 

A six-component internal force balance (AeroLab) is installed in the ADWT. The balance measurement ranges are 
±150 lb for lift, ±100 lb for side force, ±25 lb for axial force, ±250 in-lb for pitching moment, ±180 in- lb for yawing 
moment, and ±200 in-lb for rolling moment. A parallelogram linkage sting support system is able to set the angle of 
attack from -36° to 36° within the accuracy of 0.1°. A traverse system of two degrees of freedom is equipped for 
probes. 

A ceramic NACA0012 airfoil section model was used in experiments. The chord and span of the model were 10 in 
and 12 in, respectively. Along the wing span, two 1.875-in long and 6-in wide parts made of polycarbonate can be 
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mounted to the baseline model with 4 screws for clamping SETE at the trailing edge, and they can be easily 
removed to place different extended trailing edges. Several pieces of SETE with different lengths were prepared for 
tests, which were made of either aluminum or Mylar (PET-polyester). The Young’s modulus and Poisson ratio for 
aluminum were 70 GPa and 0.35, respectively, while the Young’s modulus and Poisson ratio for Mylar were 2.8 
GPa and 0.37, respectively. Figure 2 shows the NACA0012 airfoil model attached with a Mylar SETE, and the 
model parts before assembly, including the baseline NACA0012 model, two removable parts for clamping SETE, 
and typical aluminum SETE and Mylar SETE. In order to support the model, the end plates (18 by 13 in) were 
mounted to the sting balance through two aluminum bars placed in a mounting bracket. The end plates also reduced 
the three-dimensionality of flow. Figure 3 shows the model mounted on the balance for testing in the ADWT. 




(a) (b) 

Figure 2. (a) NACA0012 airfoil model with a Mylar SETE, (b) Baseline model, removable parts, and 
aluminum SETE and Mylar SETE. 



Figure 3. NACA0012 airfoil model mounted on the end plates connected to in a sting balance. 

Similarity Parameters 

The similarity parameters for a rigid airfoil with an extended trailing edge are considered. The similarity parameters 
for a rigid airfoil are Reynolds number Re c - U jz / v and Mach number M ^ , where is the freestream velocity 

and c is the airfoil chord. Since an extended trailing edge is attached to the rigid airfoil, additional similarity 
parameters should be obtained for achieving the complete flow and structure similarity. An extended trailing edge is 
generally considered to be flexible as a thin cantilever plate deforming under aerodynamic loading. The differential 
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equation for the time-independent displacement w of a thin plate is generally expressed by L[w(x,y)J = F(x,y) 
[17], where x andy are the chordwise and normal coordinates, respectively, L is a linear differential operator, and 
F(x,y) is an external distributed force on the plate. In our case, F(x,y) = Ap is the aerodynamic pressure 

difference across the plate or aerodynamic load. For a homogenous plate, the operator L - D E V 4 is the bi-harmonic 
operator, where D E - Eh 3 / 12(1 -v 2 ) is the plate flexural rigidity. The rigidity D E depends on the plate thickness 
( h ), Young’s modulus ( E ), and Poisson ratio (v). 

Let us introduce the following the non-dimensional variables F = F /q^, w = w/ 1 , x = x// , y - y // , where 
#oo = PJUt / 2 is the dynamical pressure and / is the length of SETE. Thus, we have the non-dimensional equation 
GjV 4 w = F , where the main similarity parameter for aeroelasticity is G 7 = D E /(^q^) . The functional relations for 
the lift and drag coefficients are C L = fj( a, Re c ,M ao ,G 1 , 8,1 / c) and C D = f 2 ( a, Re c ,M ao , G lf 8, l / c ) . For SETE, 
the similarity parameter G { should be sufficiently large such that the shape of SETE is not changed significantly by 
flow. For the present low-speed experiments, M ^ is very small and Re c is fixed at 4.74 x 1 0 5 . Although the 
effects of Re c and M ^ are worthwhile for further investigation, we mainly consider here the aerodynamic behavior 
of the NACA0012 airfoil model with SETE in the parametric space ( a, 8,1 / c) . 

Results and Discussions 

Static Extended Trailing Edge 


Figure 4 shows the lift and drag coefficients as a function of angle of attack (AoA) for the baseline NACA0012 
airfoil model at Re c = 4.74xl0 5 . According to McCormick’s formula [18], the lift slope depends on the aspect 
ratio (AR) by 

a _ dC L a 0 AR 

da AR + 2(AR + 4)/(AR + 2)’ K 

where a 0 - 2n according to thin airfoil theory. As shown in Fig. 4, a fit to the lift data of the baseline model using 
Eq. (1) gives the effective aspect ratio AR - 4.42 that is larger than the physical aspect ratio of 1.2. The end plates 
substantially reduced the three-dimensionality of flow. As shown in Fig. 5, the drag coefficient is represented as a 
function of the lift coefficient squared and fitted using the classical relation C D = C D0 +C 2 L 7 ne AR . Therefore, we 
obtain the zero-lift drag coefficient C D0 - 0.014 and the Oswald span efficiency e = 0.8 for AR = 4.42 . Although 

the finite end plates do not completely eliminate the three-dimensionality of flow, they modify the tip vortices acting 
like large winglets. 

In tests, an aluminum sheet (0.216 mm thick) and a Mylar sheet (0.254 mm thick) have been used as SETE, whose 
deflection angle was set a priori. The plate flexural rigidities of the aluminum sheet and Mylar sheet are 
D e =0.067 Nm and D E =0.0044 Nm , respectively. Figure 6 shows the lift coefficient as a function of AoA and 

drag polar for the model with SETE at Re c = 4.74 xlO 5 , where C L and C D are based on the projected wing area of 
the baseline model. The lengths of SETE are 5% and 10% of the chord of the baseline model, and the preset 
deflection angles of SETE are 0, 5, 7 and 10 degrees. According to Fig. 6, the C L distribution is shifted up and the 

lift is enhanced depending on the deflection angle and relative length of SETE, while the drag polar basically 
remains unchanged. The drag polar curves are collapsed for different deflection angles, indicating that the zero-lift 
drag and Oswald efficiency are not changed. Stall occurs a few degrees earlier particularly for larger deflections. 
For the model with zero deflection, C L is almost the same as that of the baseline model when AoA is less than 10°, 
but it is larger for higher AoA until stall. Typically, four and five repeated tests were conducted for each case, and 
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the results were averaged. The measurement uncertainties in C L are 0.001 for low AoA and 0.01 for high AoA. 
For C D , the measurement errors are 0.0002 for low AoA and 0.003 for high AoA. 

Further, the drag coefficient is plotted as a function of the lift coefficient squared, as shown in Fig. 7(a). All the data 
are basically collapsed into a single straight line, which is consistent with Fig. 6(b). Both the zero-lift drag 
coefficient and Oswald efficiency do not significantly change within the measurement accuracy for different 
deflection angles (small enough) and lengths of SETE. It is indicated in Fig. 7(a) that the increased C D results 

mainly from an increase of the induced drag C 2 L / neAR associated with the lift enhancement by SETE. It is 
expected that this drag penalty will be reduced as the aspect ratio increases. To examine the aerodynamic 
effectiveness of SETE, the lift-to-drag ratio (L/ D ) is plotted in Fig. 7(b) as a function of C, . The plots of L/ D 

against C, are approximately collapsed. One implication is that SETE is able to increase C, while the variation in 
L / D is relatively small. The significance of these findings is that SETE could improve the aerodynamic efficiency 
of cruise flight. 



-5 0 5 10 15 20 

AOA (deg) 

Figure 4. Lift and drag coefficients as a function of 
AoA for the baseline NACA0012 airfoil model at 

Re c = 4.74 x 1 0 s . 



CL 2 


Figure 5. Drag coefficient as a function of the lift 
coefficient squared for the baseline NACA0012 airfoil 

model at Re c - 4.74 xlO 5 . 
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(a) (b) 

Figure 6. (a) Lift coefficient as a function of AoA and (b) drag polar for the NACA0012 airfoil model with 
SETE at Re c = 4.74 x 10 5 . 
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(a) (b) 

Figure 7. (a) Drag coefficient as a function of the lift coefficient squared, (b) Lift-to-drag ratio and as a 
function of the lift coefficient for the NACA0012 airfoil model with SETE at Re c - 4.74 xlO 5 . 

Comparison with Gurney Flaps 

The aerodynamic characteristics of SETE are compared with those of Gurney flap, which is a short wall attached to 
the trailing edge that is perpendicular to the chord line on the pressure side of an airfoil, as illustrated in Fig. 8. Both 
SETE and Gurney flap increase the lift at the cost of increasing the drag. However, there are some important 
quantitative differences between SETE and Gurney flaps on the drag penalty. The measured results by Li et al. [7] 
for a NACA0012 rectangular wing indicate that L/ D as a function of C L for the wing with Gurney flaps (0.5-3%c 
heights) are significantly lower than that for the baseline wing except at higher C L . Unlike SETE, the curves of 
L/ D as a function of C L for Gurney flaps do not collapsed. Similar results can be found for other wings with 
Gurney flaps. The general conclusion is that Gurney flaps are not suitable for cruise flight due to the reduced L / D 
by the larger drag penalty. After fitting the data of Li et al. [7] using C D = C D0 + C 2 L / neAR , it is found that the 
zero-lift drag coefficient C D0 increases in a roughly linear fashion from 0.0065 for the baseline model to 0.021 for 
the 3%c high Gurney flap, where the effective aspect ratio is 46.5. The drag increase is caused by the extra form 
drag associated with open flow separation behind a Gurney flap. 

To provide a direct comparison of SETE with Gurney flaps, Gurney flap measurements were conducted using the 
same baseline NACA0012 airfoil model with the same chord Reynolds number. The Gurney flap height ranges 
from 1.2%c to 6.75%c. Figure 8(b) shows the NACA0012 airfoil model with a Gurney flap. Figure 9 shows the lift 

coefficient as a function of AoA and drag polar for the model with Gurney flaps at Re c - 4.74 xlO 5 . The C L curve 
is shifted up depending on the Gurney flap height, but C D increases considerably except for the cases with 1 .2%c 

and 2.0%c. The drag coefficient is plotted as a function of the lift coefficient squared, as shown in Fig. 10(a). The 
plots for the heights 1.2%c and 2%c are close to that for the baseline model, while other curves are considerably 
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shifted and the slope is also changed. This indicates that the zero-lift drag coefficient increases and Oswald 
efficiency varies for Gurney flaps. These results are similar to the previously published data for a NACA0012 
airfoil [7]. The lift-to-drag ratio is plotted in Fig. 10(b) as a function of C L . Unlike SETE, the L/ D curves are not 

collapsed, and they are lower than the baseline data. Figure 11 shows a direct comparison between a 10%c SETE 
with a 5° deflection angle and a 1.2%c Gurney flap. The vertical displacement of the tip of the SETE is about 
0.87%c that is close to the Gurney flap height 1.2%c. The lift enhancements in the both cases are almost the same, 
but the Gurney flap suffers from a larger drag penalty, giving a smaller L/ D than that for the SETE for all positive 

Q. 


To evaluate the net benefit of flow control in terms of the power required for cruise flight, Liu [19] has proposed the 
following weight criterion 
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( 2 ) 


where AW /W is the ratio between the control system weight AW and the total aircraft weight W , and 
rj FCr =P FCr / AW is the actuating power density. For an ideal, weightless passive control device, the criterion is 


reduced to 


6 AC " + l AC r 


7 C n 7 C T 


= g>0 . 


The function g is interpreted as a benefit margin for flow control. 


( 3 ) 


When g is positive, less engine power is 


required for cruise flight with a passive flow control device like SETE and Gurney flap, and therefore the net benefit 
is achieved. Figure 12 shows the benefit margin for SETE and Gurney flaps. For the SETE of the 10%c length and 
deflection angle of 5°, its benefit margin remains positive and approaches to zero as AoA reaches to 13°. This 
further confirms the feasibility of this SETE for improving the cruise flight efficiency. In contrast, for the 10%c 
SETE with a deflection of 14°, the benefit margin becomes negative at about 5° of AoA. The Gurney flap with the 
1.2%c height has a similar behavior, while for the Gurney flap with the 4.12%c height g is considerably below zero 
after about 3° of AoA. 



7 



3 rd International Symposium on Integrating CFD and Experiments in Aerodynamics 
20-21 June 2007 

U.S. Air Force Academy, CO, USA 




(a) (b) 

Figure 9. (a) Lift coefficient and (b) drag coefficient as a function of AoA for the NACA0012 airfoil model 
with Gurney flaps at Re c -4.74x 10 5 . 
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Figure 10. (a) Drag coefficient as a function of the lift coefficient squared, (b) Lift-to-drag ratio as a function 
of the lift coefficient for the NACA0012 airfoil model with Gurney flaps at Re c = 4.74 xlO 5 . 
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(a) (b) 

Figure 11. Comparison between SETE and Gurney flap on (a) the lift coefficient as a function of AoA, and (b) 
the lift-to-drag ratio as a function of the lift coefficient squared for the NACA0012 airfoil model at 

Re c = 4.74 x 1 0 s . 



Figure 12. Benefit margins of SETE and Gurney flaps for cruise flight. 


Thin Airfoil Theory 


Like conventional flaps, the lift enhancement by SETE is generated by a change of the camber at the trailing edge. 
For geometrical simplicity, the main wing and SETE are modeled by straight segments, as shown in Fig. 1(b). The 
total projected chord is c'=c + lcosS . The slopes of the chord line segments are 0 for 0 < x<c and -tanS for 
c < x < c ' . According to the classical thin airfoil theory assuming that the slope is small [20], the strength of the 
vortex sheet is 


f 


Y(0) = 2U w 


v 


1 + cos 6 

h 

sinO 



\ 

sin n 6 , 

/ 


(4) 


where the angular variable is related the main chordwise coordinate by x-(d / 2 )(1- cos 6 ) . The coefficients in 


Eq. (4) are A 0 = a + tan 8 ( 1 - 6 C / n ) and A n -(2/nn )tan 8 sin( nO c ) (n = 1,2,---). Here the angular variable 9 C , 
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which corresponds to the breaking point between the main wing and SETE, is given by 6 C - cos ] \(z - l)/(s + 1)\, 
where s = (l / c)cos 8 . 

The sectional lift coefficient is given by 

C { = 7 — = 2 7i( 1 + s)[a + tan 8 (l - 6 C / n + sin 0 C / tz ) ] . ( 5 ) 

(1/ 2)pJJ ao c 

The positive deflection angle of SETE shifts the C { curve upward (lift enhancement) and slightly increases the lift 
slope. The moment coefficient around the 1/4 total chord is 

M' 

e m , cW4 = , ■ . cV ; f2 2 = (1 / 4 )(1 + e) 2 tan S (sin 29 c - 2 sin 9 C ) , (6) 

where C mcV4 is defined to be positive when it makes the leading edge up. Clearly, the positive deflection angle 

produces a negative moment. Furthermore, the pressure coefficient difference between the pressure and upper 
surfaces of the airfoil is AC p =2/(x)/U ao . 


For a finite rectangular wing, the effective angle of attack is a -(1 + t )C l / n AR , where AR is the effective aspect 
ratio and the parameter r is a function of AR/ 2 n for a rectangular wing [21]. A convenient regression formula is 
r = -0.047 6( AR/ 2n) 2 +0.2195( AR/ 2tz) for AR/ 2n e [0.25,1.75 ] . Hence, the lift enhancement 


AC l -C l - C Lbaseline for a finite rectangular wing can be calculated, where the lift coefficient for the wing with 
SETE and baseline wing are, respectively, 


Q = 


7 + 


2(1 + t) 
AR 


C l9 


and 


C T 


7 + 


2(1 + t) 
AR 


2na . 


According to the above relations, the lift enhancement by SETE is 


*C l =2 


7 + 


2(1 + t) 
AR 


7i£Cc + (l + £) tan d{n -0 C + sin 0 C )\. 


(7) 


( 8 ) 


(9) 


The difference of the pressure coefficient between the pressure and upper surfaces of the main wing is sensitively 
affected by SETE. These results indicate a direct correlation of the aerodynamic loading between the main airfoil 
and SETE. Figure 13 is a comparison between the predicted AC L by thin airfoil theory and experimental and CFD 


data for the NACA0012 airfoil section (AR = 4.42 , l / c = 0.1) with SETE. The near-linear relation given by thin 
airfoil theory reasonably describes the initial trend of the experimental results for small deflections. However, the 
experimental data saturate as the deflection angle increases further, which may be caused by a reduced effective 
deflection due to viscous effect. Two-dimensional CFD calculations using a Navier-Stokes code (Fluent), assuming 
that the boundary layer on the NACA0012 airfoil with SETE is fully turbulent, predict the smaller lift enhancement. 


Thin airfoil theory has been expanded for the lift problem of an airfoil with a Gurney flap [22]. A solution of the 
generalized thin-airfoil integral equation for the vortex strength can be obtained by the method of successive 
approximations. The Gurney flap lift enhancement is interpreted as a special camber effect. The lift coefficient 
increment (pitching moment as well) is given as a square-root function of the relative Gurney flap height, i.e., 

AC l oc 4hTc . For comparison with Gurney flap, the vertical displacement of SETE is defined as 
h/c = (l/c)sinS. Therefore, AC L by SETE depends on the relative length l/c and deflection angle 5 . Figure 
14 shows the normalized AC L as a function of h/ c for SETE and Gurney flap according to thin airfoil theory. 
Interestingly, AC L increases in a linear fashion as an increase of the deflection angle for SETE with a fixed length, 
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while AC l follows a square-root relation for both Gurney flap and SETE with a variable length for a fixed 
deflection angle. 



Figure 13. Measured lift enhancement compared with 
thin airfoil theory and CFD data for the NACA0012 
airfoil model with SETE. 



h/c 


Figure 14. Comparison of normalized lift 
enhancement between SETE and Gurney flap based 
on thin airfoil theory. 


Deformation and Bending Moment 


The aluminum and Mylar sheets as SETE are usually set at nominal deflection angles in experiments assuming they 
are sufficiently rigid. Nevertheless, small aeroelastic deformation relative to the neutral position always exists. To 
estimate the elastic deformation of SETE, we consider a 2D clamped thin plate deformed due to aerodynamic load. 
The non-dimensional deformation w = w/ 1 is given by 

w(x ) = GJ 1 f dxj f dx 2 j* dx 3 f AC p (x 4 )dx 4 + Ax 3 / 6 + Bx 2 / 2 , (10) 

Jo Jo Jo Jo 

where x = x / 1 is the normalized chordwise coordinate from the breaking point between the main element and 
SETE, G 1 =D E /(l 3 qJ , A = 2G; 1 (g 1 -g 2 ) , A = G?(g 2 - 2g J , and 

r 1 r*i rx 2 _ 

gj = I dxj I dx 2 1 AC ( x 3 )dx 3 , 

Jo Jo Jo 

C 1 C * 1 _ 

g 2 - J dx 3 J AC p ( x 2 )dx 2 . 


The pressure coefficient difference AC p = 2y(x) /U^ from thin airfoil theory is substituted into Eq. (10) and 
numerical integration is carried out. Figure 15 shows the normalized deformation w/sinS as a function of 
x = x/ 1 for the 0.216 mm thick aluminum sheet (G 1 =5.74 , 10%c) and 0.254 mm thick Mylar sheet (G 1 =0.38 , 
10%c). The deformation of the aluminum SETE is less than 1% of the vertical tip displacement l sin 8 , while for 
the Mylar SETE it could be 13%. To estimate the actuating power required for deflecting SETE, the bending 
moment relative to the breaking point between the main element and SETE is estimated. The bending moment 
coefficient is defined as 


C 


M b 



AC p (x)xdx . 
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Figure 16 shows the bending moment coefficient as a function of AoA for SETE with a deflection angle of 5 
degrees. The bending moment for the 2D case is in the order of qj 2 . 



Figure 15. Normalized deformation of the aluminum SETE and Mylar SETE with a deflection angle of 5 
degrees for different AoAs. 



Figure 16. Bending moment of SETE with a deflection angle of 5 degrees relative the breaking point between 
the main element and SETE. 

CFD SunuMions 


In order to provide further insights to the lift enhancement by SETE, the flow characteristics around a NACA0012 
airfoil with SETE were studied using CFD simulations Re c = 4.74 xlO 5 . The problem of a steady turbulent flow 

over the airfoil was solved using the commercial code Fluent on a C-grid. Grid independence was validated by 
comparing pressure distributions along the airfoil at two different angles of attack with published results [23]. The 
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lift and drag coefficients were also compared with the published data [23]. The boundary layers were assumed to be 
fully turbulent, and the realizable k-epsilon turbulence model with the non-equilibrium wall function was used. 

Figure 17 shows how the defection angle of SETE affects the pressure coefficient distribution on the airfoil for zero 
AoA. The aerodynamic loading on SETE itself (x/c>7)is not large, and the lift enhancement is mainly produced 
by the overall aerodynamic change on the main airfoil induced by SETE. Comparisons of the calculated and 
measured lift coefficients are shown in Fig. 18 as a function of AoA. The calculated results for the 2D airfoil are 
converted to the 3D wing results using Eq. (7). The calculated and experimental results are consistent until the 
calculated lift coefficients level at higher AoAs. As indicated in Fig. 13, the calculated increment in the lift 
coefficient exhibits the same trend as the experimental and theoretical results as a function of the deflection angle of 
SETE. Figure 19 shows velocity magnitude contours around the NACA0012 airfoil with zero AoA for 0° and 5° 
deflection angles. The asymmetry of the flow field is induced by SETE and the wake is turned downward, 
indicating a deflected momentum stream tube and the generation of additional lift. The wake structure is not 
appreciably altered although it is deflected; this implies that the parasite drag is not significantly affected. 

Further, a 10%c SETE with a 5° deflection angle is compared with a conventional plain flap that is deflected at the 
90% chord of a NACA0012 airfoil. Figure 20 shows the computed lift and drag coefficients as a function of the 
deflection angle for the SETE and corresponding plain flap. Compared with the plain flap, the lift enhancement by 
the SETE is larger while the drag coefficient is smaller. Here the lift and drag coefficient for the SETE is defined 
based on the main airfoil chord c. Even when these coefficients are based on 1.1c for the SETE, the above results 
are still qualitatively valid. The pressure coefficient distributions on the main airfoil for the SETE (5° deflection), 
plain flap (5° deflection) and Gurney flap are shown in Fig. 21, indicating that the larger pressure difference on the 

main airfoil is induced by the SETE. The Gurney flap height is 0.1csin( 5° ) « 0.87%c , which is the same as the 

maximum vertical displacement of the SETE. Figures 22 and 23 show the lift coefficient as a function of AoA and 
the drag polar for the SETE, plain flap and Gurney flap. Among these trailing edge devices, the SETE generates the 
largest lift increase at a smallest drag penalty. It should be emphasized here that the absolute accuracy of the drag 
predicted by CFD is not critically examined although the computed drag coefficient for the NACA0012 airfoil is 
close to C D =0.02- 0. 03 given by measurements for the corresponding airfoil section. Therefore, the above 

comparison in the drag polar is valid only in a relative sense. Furthermore, for the 2D airfoil, the drag polar is not 
parabolic since no induced drag is included. 


1 1 1 

O Deflection. 0 deg 

+ Deflection. 5 deg 




-0.4 ^ 1 J 1 1 L 

0 0.2 0.4 0.0 0.0 1 


x/c 

Figure 17. Calculated pressure coefficient 
distributions on the NACA0012 airfoil for the SETE 
deflection angles of 0 and 5 degrees at zero AoA. 



Figure 18. Comparisons of the calculated and 
measured lift coefficients as a function of AoA for the 
deflection angles of 0 and 5 degrees. 
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(a) (b) 

Figure 19. Velocity magnitude fields of the NACA0012 airfoil at zero AoA, (a) zero deflection of SETE, (b) 5 
degrees deflection of SETE. 




Figure 20. Comparison of the computed lift and drag coefficients between SETE and conventional flap for a 
NACA0012 airfoil at zero AoA. 



Figure 21. Comparison of the computed pressure coefficient between SETE, conventional flap and Gurney 
flap for a NACA0012 airfoil at zero AoA. 
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Figure 22. The computed lift coefficient as a function Figure 23. The computed drag polar for SETE (5-deg 
AoA for SETE (5-deg deflection), Gurney flap and deflection), Gurney flap and conventional flap (5-deg 
conventional flap (5-deg deflection) for a NACA0012 deflection) for a NACA0012 airfoil, 
airfoil. 

Conclusions 

Experiments and calculations presented show that a static extended trailing edge (SETE) attached to a NACA0012 
airfoil model is able to enhance the lift while the zero-lift drag is not significantly increased. The mechanism of lift 
enhancement by SETE is the camber effect, which is the same as that for other high-lift devices like Gurney flap and 
conventional flap. However, compared with Gurney flap and conventional flap, SETE generates a larger lift 
increase at a smaller drag penalty since it is embedded in the wake of the main airfoil. Therefore, SETE has a 
promising potential for improving the efficiency of cruise flight. The benefit margin of SETE for cruise flight is 
evaluated in comparison with Gurney flap, and the feasibility of SETE for lift enhancement in cruise flight is 
demonstrated at small angles of attack and deflection angels of SETE. The mechanical simplicity of SETE allows 
direct application to aircraft without changing the basic aerodynamic configuration. Furthermore, the deformation 
and pitching moment of SETE due to aerodynamic loading are estimated, and actuating SETE is feasible since the 
pitching moment on it is small. Future research will focus on SETE actuation using smart material like 
piezoactuators and shape memory alloys for steady and unsteady flow control. 
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